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I. INTRODUCTION 

The discovery of superconductivity in CePtaSi (Ref. Q]) 
has renewed interest, both experimental and theoretical, 
in the properties of superconductors without inversion 
symmetry. The list of such superconductors has been 
steadily growingand now also includes Ulr (Ref. 2j), 
CeRhSi 3 (Ref. CeIrSi 3 (Ref. 0), Y 2 C 3 (Ref. E), 
Li 2 (Pd 1 _ ;I ,,Pt :E )3B (Ref. and other materials. 

A distinctive feature of noncentrosymmetric crystals is 
that the spin-orbit (SO) coupling qualitatively changes 
the nature of single-electron states, namely it lifts spin 
degeneracy and splits energy bands almost everywhere in 
the Brillouin zone. This has important consequences for 
superconductivity. In the limit of strong SO coupling, 
the Cooper pairing between the electrons with opposite 
momenta occurs only if they are from the same nondc- 
generate band. According to Ref. 0, this is what hap- 
pens in CePtaSi, where the SO band splitting ranges from 
500K to 2000K and thus exceeds the critical temperature 
T c = 0.75K by orders of magnitude. The same is likely 
to be the case in other materials, for instance Li 2 Pd3B 
and Li 2 Pt 3 B, see Ref. & 

The pairing interaction in the strong SO coupling 
case is most naturally introduced using the exact band 
states;^ : 10 ' 11 which take into account all the effects of 
the crystal lattice potential and the SO coupling. In the 
band representation, the superconducting order param- 
eter is represented by a set of complex functions, one 
for each band, which are coupled, e.g. due to the in- 
tcrband tunneling of the Cooper pairs or the impurity 
scattering. 12 Since the bands are nondegenerate, the pair- 
ing symmetry is peculiar: While the Pauli principle dic- 
tates that each order parameter is an odd function of 
momentum, the gap symmetry, in particular the loca- 
tion of the gap nodes, is determined by one of the even 
representations of the crystal point group. 

If the band splitting is smaller than the superconduct- 
ing critical temperature, then the effects of the SO cou- 
pling can be treated perturbatively, in particular the pair- 
ing Hamiltonian can be constructed using the basis of the 
pure spinor states unaffected by the SO coupling. This 
approach to the theory of noncentrosymmetric supercon- 
ductivity was introduced in Ref. [H and developed fur- 
ther in Refs. GJQJ. 

One of the most peculiar properties of noncentrosym- 



metric superconductors is a strongly anisotropic spin 
susceptibility with a large residual component at zero 
temperature^ i 13 ' 14 i 16 ' 17 i 18 ' 19 i 20 The latter comes from the 
field-induced virtual interband transitions and is almost 
unchanged across T c . Its magnitude depends on the 
electron structure and is considerably smaller than the 
normal-state susceptibility, see Sec. Ill Al below. On the 
other hand, the experiments in CePtaSi (Ref. |21| ) have 
detected essentially no change in the Knight shift below 
T c , for all directions of the applied field. One possible 
explanation is that the relative magnitude of the inter- 
band contribution increases due to the combined effect of 
interactions and peculiarities of the electron structured 

The goal of this article is to explore an alternative 
mechanism whereby the residual spin susceptibility is en- 
hanced by disorder. It is well known that the SO impu- 
rities play an important role in usual, i.e. centrosymmet- 
ric singlet, superconductors, in which experiments also 
revealed a large residual Knight shift, see, e.g. Ref. I22I 
and the references therein. This is incompatible with the 
predictions of the Bardeen-Cooper-Schrieffer (BCS) the- 
ory, even when the SO coupling with the crystal lattice 
potential is taken into account^ and can be explained 
by the presence of spin-reversing scattering at the sam- 
ple boundaries or impurities i 24 ' 25 ! 26 We calculate the spin 
susceptibility of a noncentrosymmetric superconductor in 
the presence of scalar disorder. In contrast to Refs. [l3l 
and I la. we focus on the strong SO coupling limit and em- 
ploy the band representation of the pairing Hamiltonian. 

The article is organized as follows: In Sec. [IT] we calcu- 
late the spin susceptibility in the clean case. In Sec. IIIIl 
we introduce the impurity scattering in the band repre- 
sentation and derive general expressions for the suscepti- 
bility in the disordered case. In Sec. IIV1 we calculate the 
residual susceptibility in some simple models, in which 
the effects of impurities can be worked out analytically. 
We have also included an Appendix with the derivation 
of the effective SO Hamiltonian for band electrons. 



II. CLEAN CASE 

In a noncentrosymmetric crystal with SO coupling the 
electron bands are nondegenerate. The formal reason is 
that without the inversion operation one cannot, in gen- 
eral, have two orthogonal degenerate Bloch states at the 
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same wave vector k. In the limit of zero SO coupling 
there is an additional symmetry in the system - the in- 
variance with respect to arbitrary rotations in spin space 
- which preserves two- fold degeneracy of the bands. Let 
us consider a single band of electrons with the disper- 
sion given by e(k), and turn on the SO coupling with the 
crystal lattice. The system Hamiltonian can be written 
in the following form: 



Hn = 



^2 [e(k)S a p + 7(fc)er Q/3 ]aJ CQ a fe/3 . 



(1) 



k,a/3 



Here a, (3 =T,I label the spin states, a are the Pauli 
matrices, and the sum over k is restricted to the first 
Brillouin zone. The "bare" band dispersion satisfies 
e(— k) = e(k), e(g~ 1 k) = e(k), where g is any operation 
from the point group G of the crystal. The SO coupling 
is described by the pseudovector function ~f(k), which 
has the following symmetry properties: -f(k) — —j(—k), 
5 , 7(<7 _1 fc) = 7(fc). The derivation of Eq. ([I]) is outlined 
in the Appendix. 

For the tetragonal group G = C4„, which describes 
the point symmetry of CePt 3 Si, CeRhSi3 and CelrSia, 
the SO coupling has the following form: 



7(fe) =7±[0E,«(fe) x z]+i\\<f>A a ,v{k)z. 



(2) 



Here 7j_ and 7|| are constants, and (f>E,u and 4>a 2 ,u are 
the odd basis functions of the irreducible representations 
E (two-dimensional) and (one-dimensional), respec- 
tively. The representative polynomial expressions for the 
basis functions are 4>E,u(k) oc (k x ,k y ) and <f>A 2 ,u(k) oc 
k x k y k z (k x — ky). A particular two-dimensional case of 
Eq. with 7(fe) = J±(k x z) is known as the Rashba 
model— and has been extensively used to study the ef- 
fects of SO coupling in semiconductor heterostructures. 
For the cubic group & — O, which describes the point 
symmetry of L^Pdi-jcPtr^B, the SO coupling has the 
form 



7(fe) = 7o0F 1 ,«(fe), 



(3) 



where <t>F lt u are the odd basis functions of the vector 
representation Fx, e.g. 4>p liU (k) cx (k x , k y , k z ). 

The Hamiltonian fl} is diagonalized by a unitary trans- 
formation: 



a-ka = u a \{k)c k \, 



(4) 



where 




(5) 



and 9\ are arbitrary (in general fc-dependent) phases. 
The free-electron Hamiltonian then becomes 



(6) 



where £\(k) = e(fe) + A|7(fe)| = k) describes the 
quasiparticle dispersion in the Ath band. Thus the SO 
coupling lifts the spin degeneracy of the electron bands, 
so that the Fermi surface consists of two pieces defined 
by the equations £±(fc) — 0. The magnitude of the band 
splitting is given by 2|7(fe)| and might vanish, for sym- 
metry reasons, along some directions or at some isolated 
points in the Brillouin zone. We assume that the band 
structure is such that the zeros of -f(k) are not located on 
the Fermi surface. This is the case in the model ([3]), and 
also in the model if the Fermi surface is a cylinder 
around the z-axis. 

The Zeeman coupling of the electron spins with an 
external magnetic field is described by 



H z = -VbB o-a/jaL 



fe,Q/3 



= -B ^2 m \x'{k)c { kX Ck\', 



(7) 



fc.AA' 



where /is is the Bohr magneton. The components of 
rh(k) = figiv (k)&u(k) have the following form: 



m x = fi B 



lx -e 

7-L 



~ _ / ly ~ e 

m v ~ t iB I i§ 2y%+Hz 
\ 7x 



-id lulz+hy 

7-L 

-i§ 7y7^~'7m 



1± 



(8) 



m z = hb 



where 7 = 7/I7I, 7± = +7^, and 6 = 9+ - 9-. 

The expectation value of the spin magnetic moment of 
an electron from the Ath band with the wave vector k 
is m\(k) = A^s7(fc). We shall see that, although the 
interband matrix elements of rh(k) contain the arbitrary 
phases of the Bloch spinors, those will not affect observ- 
able quantities. 

Finally, we introduce the pairing interaction between 
electrons in the Cooper channel, using the basis of the 
exact eigenstates of the noninteracting problem: 



X C-k'+q/2,\'Ck'+q/2,\', (9) 



kk' q XX' 



k X=± 



where V is the system volume. The pairing potential 
satisfies the relations V\\'(—k,k') — —V\\>(k,k') — 
V\\/(k, —k'), which follow from the anti-commutation of 
the Fermi operators. The diagonal elements of the ma- 
trix V describe the intraband Cooper pairing, while the 
off-diagonal ones correspond to the pair scattering from 
one band to the other. We assume, in the spirit of the 
BCS theory, that the pairing interaction is nonzero only 
inside the thin shells of width e c in the vicinity of the 
Fermi surfaces, i.e. when |£x(fe)|, |£a'(^')I < £ c- Wc 
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further assume that it can be represented in a factor- 
ized form: V xx >(k,k') = t x {k)t x ,(k')V X y(k,k'), where 
t\(k) = —t\(—k) are nontrivial phase factors which ap- 
pear in the expression for the time-reversal operation for 
nondegenerate bands: K\k\) = t\{k)\ — k, A)^^ and 



dr 



V xx ,(k,k') = -V xx ,J2MkWy, a (k>) (10) 



0=1 



is invariant under the point group operations: 
V xx <(g~ 1 k,g- 1 k') = V xx ,(k,k'). The coupling con- 
stants V XX i form a symmetric positive-definite 2x2 
matrix, and <j) X ^ a (k) are even basis functions of an 
irreducible dr-dimensional representation Y of G. 28 
While +iQ (fe) and (/>_ jQ (fc) have the same symmetry, 
their momentum dependence does not have to be exactly 
the same. The basis functions are nonzero only inside 



the BCS shells and are normalized: 



= (>ab, 



where the angular brackets denote the averaging over 
the Fermi surface in the Ath band. 

Treating the pairing interaction ([9]) in the mean-field 
approximation, one introduces the superconducting or- 
der parameters, which in the uniform case considered 
here depend only on k and have the form A x (k) — 
t x (k)A x (k), where A x (k) = A x (—k) can be represented 
as follows: 



(11) 



Thus our system is formally equivalent to a two-band 
superconductor, in which the order parameter has 2dr 
components given by the expansion coefficients Tj X a . The 
order parameter components depend on temperature and 
vanish in the normal state. 



A. Spin susceptibility 

Writing the Zeeman Hamiltonian ([7]) in the form Hz = 
—A4B, where M. is the operator of the total spin mag- 
netic moment of electrons, we define the magnetization 
as M = V~ 1 (M) (the angular brackets denote the ther- 
modynamic average). In a weak field, Mi — Y^j XijBj> 
where Xij is the spin susceptibility tensor. 

Introducing the four-component Nambu operators in 
the band representation: Ck = (ck+, Qc-, c_ k , , c_ k + ) T ', 
we combine the normal and anomalous Green's 
functions 29 into a 4 x 4 matrix Green's function 



g(k 1 ,k 2 ;r)=-(T T C kl (r)CUo)). 



(12) 



In the clean system, the Green's function is diagonal in 
momentum: 



G(k,w n ) = 



G(k,LO n ) 

-Ft(fc, w „) 



-F(k,cj n ) 
-G T (-k,-uj, 



(13) 



Here ui n — (2n + 1)ttT is the fermionic Matsubara fre- 
quency (in our units kg — 1), and the hats here de- 
note 2x2 matrices in the band space. In the thermody- 
namic limit V — > 00 the magnetization can be expressed 
in terms of the Green's functions as follows: 



/d 3 k 
^-^tY7h(k)G(k,U n ), 



(14) 



where the matrices rh are given by Eqs. {5J. 

In the presence of magnetic field we have G{k, u> n ) — 
[G ' 1 (k,u! n ) - Y.z(k)]^ 1 , where 



G 1 {k,uj n ) = 



u n - i(k) -A(fe) 
-At (As) iw„ + £(k) 



(15) 



A XX /(k) = 5 XX 't x (k)A x (k), and the Zeeman coupling is 
described by 



^ z{k) ~ { rh T {-k)B 
Expanding Q in powers of B, we obtain: 
d 3 k 



(16) 



Xij 



(2tt) s 



tiijhiGrhjG — rhiFrhjF'j , (17) 



where m(k) — m T (~k). The Green's functions here are 
calculated at zero field and have the following form: 

G XX '(k,uj n ) = 5xx'G x (k,u n ), 

F X y(k,UJ n ) = 8 X yt X (k)F X {k,UJ n ), (18) 

Fl x ,(k,u n ) = 5xx>tl(k)FZ(k,u n ), 



where 



G\(k, uj n ) 
Fx(k,Lu n ) ■ 



^ 2 +£ 2 (fc) + |A A (fc)| 2 ' 



Ax(k) 



(19) 



c 2 + e 2 (fc) + |A A (fe)i 2 



Inserting the expressions ||5J) for the electron magnetic 
moment in Eq. (| 1 T[) and using the identities 



mi,\\> = -t x t x >m iiX y, 
m h++ m h++ = m h — m h — 
m iy+ _m j _ + = n 2 B (S i: j - 7^ + ie ijk %), (20) 



we find that the susceptibility tensor can be represented 
in the following form: 



ij + Xij ■ 



(21) 



Here 



Xij 



-Mb 



707*7,(G 2 a + |A| 2 ) 



ANx {liljYx)) 



(22) 
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are the intraband contributions, which are determined 
by the thermally-excited quasiparticles near the Fermi 
surfaces, N\ is the density of states in the Ath band, and 



cosh 2 ( 



|A A (fe)|V2T) 



(23) 



is the angle-resolved Yosida function. 
The last term in Eq. §TI$) is 



Xij 



(G+G- + ReF+F_) 



(24) 



In the normal state this becomes 



"Mb 



d 3 k S, 



(2tt)3 | 7 | 



[/(£+)-/(£-)], (25) 



where /(e) = (e e / T + 1) _1 is the Fermi-Dirac distribution 
function. In contrast to the intraband susceptibilities 
X^j, which depend only on the quasiparticle properties 
in the vicinity of the Fermi surfaces, Xij is determined 
by all quasiparticles in the momentum-space shell "sand- 
wiched" between the Fermi surfaces. We therefore call 
this the interband contribution (its physical origin is dis- 
cussed in Sec. IIIBI belowh A straightforward calculation 
shows that even at T = 



Xij \A^0 Xij\A= 
Xij\ A=0 



l7l : 



inm«i, 

A 



which means that Xij is almost unchanged when the sys- 
tem undergoes a phase transition in which only the elec- 
trons near the Fermi surface are affected. 

Thus we arrive at the following expression for the spin 
susceptibility of a clean superconductor: 



Xij = Xij + 



(26) 



(N + +N_)/2, and 
= 1 and 



where Xij is given by Eq. (|25| . Np = 
px = Nx/Np. In the normal state, Y\ 

XN,tj = Xij + p%N F ^2 px (itjj)x ■ (27) 

A 

At zero temperature, there is no excitations (Yx = 0) and 
the intraband contributions are absent, but the suscepti- 
bility still attains a nonzero value given by Xij ■ The tem- 
perature dependence of the susceptibility in the super- 
conducting state at < T < T c is almost entirely deter- 
mined by the intraband terms, with the low-temperature 
behavior depending crucially on the magnitude of the SO 
coupling at the gap nodes. 19 While in the fully gapped 
case the intraband susceptibility is exponentially small 
in all directions, in the presence of the lines of nodes it is 
proportional to either T or T 3 , depending on whether or 



H 

X 



0.8 



T/T c 



FIG. 1: The temperature dependence of the transverse com- 
ponents of the susceptibility for the 2D model (the solid line), 
and of all three components for the 3D model (the dashed 
line). The \ Z z component in the 2D case is temperature- 
independent. 



not the zeros of A\(k) coincide with those of 7(fc), see 
Ref. [H for details. 

Further steps depend on the pairing symmetry and the 
electron structure of the superconductor. Let us evaluate 
the expression (f26|) in the case when the SO coupling 
is small compared to the Fermi energy, i.e. | — y | <C ep, 
and the gaps in both bands are isotropic and have the 
same magnitude. It is legitimate to neglect the difference 
between the densities of states in the two bands: p + = 
p- = 1, and we also have Y+(k,T) = T_(fc,T) = Y(T). 
The interband contribution (|2"ST) is reduced to 



Xij = 2p 2 B Np(5 l:j - (jiJj) F ), 



(28) 



where ((•••)) f denotes the average over the Fermi sur- 
face defined by the equation e(k) = 0. The normal- 
state susceptibility is isotropic: XN,ij — XP&iji where 
Xp = 2p B Np is the Pauli susceptibility, while in the su- 
perconducting state 



Xtj = Xij + 2HB N F(%lj)FY(T). 



(29) 



Suppose the Fermi surface is a cylinder around the 
z axis and the SO coupling is described by the Rashba 
expression: j(k) = ^±(kx z). We will refer to this as the 
two-dimensional (2D) model. The susceptibility tensor is 
diagonal, with 



1, 



X X *=Xyy = ^ + Y(T)} X P, 



Xp- 



(30) 



For a spherical Fermi surface in a cubic crystal with 
the SO coupling described by ~f(k) — jok [the three- 
dimensional (3D) model], the susceptibility tensor is 
given by 



X.xx — Xyy Xzz 



Y(T)} XP . 



(31) 
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The temperature dependence of Xij f° r the 2D and 3D 
models is sketched in Fig. [TJ In both cases the suscepti- 
bility at T = is nonzero but still considerably less than 
in the normal state. 



B. Origin of the residual susceptibility 

The mean-field pairing Hamiltonian in 
the band representation is given by H sc = 

(l/2)E fc ,A[AA(fe)4 A c t _ fc ,A + H.c.]. This can be 
transformed into the spin representation using Eqs. (f?|): 

Hsc = (1/2) Eh,afi[ A ^i k ) a L a -h,p + H - c -L where 

A„/3(fc) = ^M QA (fc)A A (fc)u/3A(-fc) 
A 

= -[u{k)k(k)u\k){ia 2 )] al3 . (32) 

To obtain this, we used the identity Ufj\(—k) = 
X^* ^) f3y u y\(k) , which follows from the definition 

of t\(k). The next step is to write A = (A + + A_)f /2 + 
(A + — A_)r3/2, where f are the Pauli matrices. Insert- 
ing this in Eq. (|3"2"j) and using 

u(k)r 3 u\k) =j(k)T, (33) 

we obtain the order parameter in the spin representation 
as follows: 

A Q/3 (fc) = tp(k)(ia 2 ) a p + d(k)(iaa 2 ) a p, (34) 

where 

m = - A+W ; A - W (35) 

is the spin-singlet component, and 

i{h) _ _ A +W -A-( fc ) , m (36) 

is the spin-triplet component. 30 In agreement with the 
results of Ref. Q, only one component of the triplet 
order parameter survives (is "protected" ) in the limit of 
a large SO band splitting. 

We see that the residual susceptibility is not related 
to the presence of the spin-triplet order parameter: Ex- 
pression (I36|) vanishes if the both gaps are the same, but 
the residual susceptibility is still nonzero and given by 
Eq. (|25p • The origin of the residual susceptibility can be 
understood using the following argument.— For a given 
wave vector k, there are two electron states |fc, ±) with 
energies £,±(k), in which the expectation values of the 
spin magnetic moment are m±(k) — ±fig / y(k). Due to 
the presence of the projector 8^ — %*fj, only the compo- 
nent of B which is perpendicular to m±(k) contributes 
to the susceptibility (|23|) . Calculating the energies of the 



states \k, ±) in the second order of the perturbation the- 
ory with the help of the identity (f2"TJ)) . we obtain: 

To find the total energy shift, we sum the contributions 
from all k and both bands, with the weights given by 
the Fermi-Dirac distribution. In this way we recover Eq. 
(|25[) . Thus the intcrband susceptibility appears, similarly 
to the Van Vleck paramagnetism of atoms and solids, 
in the second order of the perturbation theory, due to 
the field-induced virtual transitions between the SO split 
electron bands. 



III. EFFECTS OF IMPURITIES 

The effect of scalar impurities on the spin susceptibility 
in the superconducting state is described by the Hamil- 
tonian H = Ho + Hz + Hi n t + Hi mp , where the first three 
terms are given by Eqs. ©, ([7]) and ©, and 

H imp = I d 3 rJ2Uim P (r)^ a (r)^ a (r). (37) 

The disorder potential Ui mp (r) is assumed to be a 
random function with zero mean and the correlator 
(U imp (ri)Ui mp (r 2 )) = n imp U 2 S(r 1 - r 2 ), where 
the impurity concentration, and Uq is the strength of an 
individual point-like impurity. The electron field opera- 
tors have the form 

ip a {r) = -^=^2u aX {k)e tkr c kx , 

from which we obtain the band representation of the im- 
purity Hamiltonian: 

H imp = ^ ^ ^ Uimpjk - k')w\\> (fe, fc')4 A c fe /A' . (38) 

fcfc' AA' 

Here Ui mp {q) is the Fourier transform of the impu- 
rity potential, (U imp (q)U imp (q')) = n imp U^VS q - q ' , and 
w(k,k') — u> {k)u{k') — w'(k',k). We see that the im- 
purity scattering amplitude in the band representation 
acquires both intraband and interband contributions and 
also becomes anisotropic, even for isotropic impurities. 12 

Averaging with respect to the impurity positions re- 
stores translational invariance: (Q(ki, k 2 ', ^ n ))imp = 
<5fei,fe 2 ^(' c ' w «)- The disorder-averaged Green's function 
here has the same matrix structure as in the clean case, 
see Eq. (1131) . The magnetization in the superconducting 
state is determined by Eq. (fT4|) . in which the Green's 
function of the clean normal metal should be replaced by 
its disorder average. 

The average matrix Green's function satisfies the 
Gor'kov equations: (Gq 1 — ^imp — ^z)Q — 1, where Go 
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is the Green's function at zero field in the absence of im- 
purities, given by Eq. (Ti~5"]) . the impurity self-energy in 
the self-consistent Born approximation is 



where 



f d 3 k' 1 . ^ r d 3 k ■ 

J j^W(k,k')g(k',u n )W(k',k), (39) ^W-^oE] (2^ FA(fc ' 



(43) 



(27T) 

and is the Zeeman self-energy (TIB)) . The 4x4 matrix 
PF is defined as follows: 



It is straightforward to show that [w T (— k', — k)]w = 
tl(k)t x ,(k>)w xx >(k,k'). 



A. Zero-field solution 

In the absence of magnetic field, £z = 0, and, since 
Q is band-diagonal, we seek the solution of the Gor'kov 
equation in a band-diagonal form given by Eqs. (|18jl. For 
consistency, we require that the Nambu matrix compo- 
nents of the self-energy are also band-diagonal: 

(T imp )f x , (fe, Wn ) = 5^A' Q a x b (k)af(k, «„), (40) 

where a, 6 = 1,2 are the Nambu (particle-hole) indices, 
and gii(fc) = gf(fe) = 1, e! 2 (fc) = i A (fc), ^(fe) = 
t A (fc). From Eq. flU we obtain: 

CT n (fc,w„) <7 12 (fc,w n ) 
a 21 (k,uj n ) & 22 {k,uj n ) 

X ( U ( ^ fc) . (41) 

Using the identity (j3"3")) one can show that the integrands 
on the right-hand side of these equations have the follow- 
ing form: 

u{k')G{k\uj n )u\k') 

G + {k',uj n ) + G^{k',uj n ) . 
= 2 T ° 

+ G ' +(fc> " ) - G - (fc> " ) 7(fc0f i (42) 

ete. Assuming that G\ and Fa are even functions of 
momentum (the self-consistency of this will be verified 
below) , the last line in Eq. (|4"!?]) is odd in k' and therefore 
vanishes after the fc'-integration. Then, 

^(fc,^) = -af 2 (k,-uj n ) = Si(w„), 
a{ 2 (k,uj n ) = cr A **(fe,o;„) = S 2 (o; n ), 



Absorbing the real part of £ into the chemical potential, 
we have £i(o; n ) = iSi(w„), where Si is odd in u> n . 

Solving the Gor'kov equations we obtain the disorder- 
averaged Green's functions: 



G x (k,u n ) 



F x (k,uj n ) 



i& n + i\(k) 



u ) 2 n + e x (k) + \D x (k,u )n )\ 2 ' 



D x (k,u n ) 



(44) 



Q 2 +e x (k) + \D(k,u Jn W 



where ui n = ui n — Ei(w„) and D x (k,u n ) — A x (k) + 
^2(^n)- Substituting these expressions into Eqs. (f43j) . 
we arrive at the self-consistency equations for the renor- 
malized Matsubara frequency and the gap functions: 



xX /C J 2 + \D x (k,u; n )\ 2 / x 



D x (k,u; n ) = A x {k) 



E 



Dx'(k,w n ) 



y/u% + \D X ,{k,W n )\* „ 



(45) 



.(46) 



Here T = l/2r is the elastic scattering rate, and r = 
(2imi mp UQN f) -1 is the electron mean free time due to 
impurities. One can now see that, while w„ is odd in u n , 
D x are even in both k and oj n , so that our assumptions 
are self-consistent. 



B. Spin susceptibility 

At B 7^ 0, both the Green's function and the impurity 
self-energy acquire field-dependent corrections: Gb^o = 



imp 



ST, 



imp • 



Here Q and £ 



imp 



are the average Green's function and the impurity self- 
energy at zero field found in Sec. IIII Al Treating the 
Zeeman coupling as a small perturbation, we find 5Q = 
G(T. Z + ST, imp )g. 

Magnetization M is given by Eq. (fT4)) . The contri- 
bution from the zero-field Green's function vanishes af- 
ter the momentum integration, meaning that there is no 
spontaneous magnetism of the Cooper pairs in the super- 
conducting state. Thus the magnetization is determined 
by 6G and can be written as M = Mi + M2 , where 



Mi = 



d 3 k 



trm(GY, z g) 
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(47) 
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corresponds diagrammatically to a "bubble" containing 
two disorder-averaged Green's functions, while 



Mo 



d 3 k 

(27T) 



3 trm(^E lmp g) n (48) 



describes the impurity vertex corrections (recall that the 
upper indices here label the Nambu matrix components) . 

The contribution to the susceptibility from Mi can 
be calculated similarly to the clean case in Sec. Ill Al 
Xi,ij = dMx,i/dBj = xfj + Xij + Xij, where x}j are the 
intraband and Xij the intcrband susceptibilities. Since 
the latter is determined by all quasiparticles between the 
"+" and "— " Fermi surfaces, its change below the super- 
conducting transition is negligibly small. One can show 
that the interband susceptibility in the normal state is 
not sensitive to impurities, therefore Xij is given by its 
clean normal state expression (|2~51) . The intraband sus- 
ceptibilities are given by 



Xij 



2 rpV^ f ^k, - - 



\FA 2 ) 



Mb^a (jilj) 



R+|D A |2)3/ 5 



, (49) 



The Green's functions in the first line here are given by 
the disorder-averaged expressions (|44|) . In contrast to 
the clean case, however, it is not possible to calculate 
the Matsubara sums before the momentum integrals. In 
order to do the momentum integrals first, one should 29 
add and subtract the normal-state intraband susceptibil- 
ity, which is not affected by impurities, see Eq. (f2"T)) . 

Let us now calculate the impurity vertex corrections 
(gSJ). Substituting 8Q = £(£ z + <5£ lmp )£ into Eq. 
we obtain the following equation for 8Y>i mp : 

f d 3 k' 

8Y, imp (k,uj n ) -riimpUo J —^W{k,k')g(k',u) n ) 
x<5£ imp (fc', u n )G(k', u> n )W(k', k) 



TT 2 I dW_ 



W{k,k')g(k',u n 



(2nf 

x£ z (fc')S(fc', uj n )W{k',k). (50) 

The contributions to iW 2 from the interband components 
of 5Tii mp can be neglected compared to those from the 
intraband components, because the former contain the 
momentum integrals of the products of the Green's func- 
tions from different bands, which are at least by a factor 
of |A|/|7| <C 1 smaller than their same-band counter- 
parts. One can seek the field-induced correction to the 
impurity self-energy in the band-diagonal form similar to 
Eq. (SOD: 

{SZ imp )f x , (k, w„) = 6 X y gt b (k)Sa a x b (k, «„). (51) 

For the same reason, one should retain only those terms 
on the right-hand side of Eq. (|5"D|) in which both 



Green's functions have the same band index. Since 
mx\(—k) = —mxx(k) = -XfiBj(k), see Eqs. ([5]), we 
have (E z )ft(fc) = (Sz)f A (fc) = -\» B j{k). Inserting 
this and the expressions (|5Tj) in Eq. ([50]) and using the 
identity 

we obtain the Nambu matrix components of the right- 
hand side of Eq. (J5T 



(R.H.S.)f(k,u n ) = -\nBQ a x b (k)J2Uk)B 3 

1 [ d 3 k' 

where i? A x = G\ + |F A | 2 , Rf = (G A - G A )F A , Rf = 

(Gx - G X )F*, Rf = G\ + \F X \ 2 , and G A (fc,^„) = 

Gx(—k,—u) n ). This suggests that one can seek Sa x b in 
the following form: 

foJH*)^) = Sal 2 (k,uj n ) = Xj(k)X(uj n ), 



5<r{ 2 (k,u; n ) = -8a x '*(k,u n ) = Xj(k)Y(u n ). 



(52) 



After some straightforward algebra, we obtain from 
Eq. ([50)1 the linear equations for X(oj n ) and Y(oj n ): 



X t - + A %ii Yi + A* %ij Y*) = X , 

j 

Y t - £(2A^ + A 3iij Yj + A^Y*) = Y , 



(53) 



The notations here are as follows: 



\Dx? 



(«a + |D A |a)3/a 



A I \ 1 I - - i^nPx 

A 2<ij (Un) = 7L^ V^(^ + |^ a | 2) 3/2 

a i \ r Y- A- 2cD 2 + |£ A | 2 
4wK) = j 1_,PX \^7i ( -, + |jDA|2)3/2/ a 

A ii%j (w n ) = - l^Px.JUjjz 



X Q .i(uj n ) = -HB/jAx^u^Bj, 

3 

Yo,i{oJ n ) = -2/xb ^ Al di (u n )Bj . 



We see that, while X{uo n ) is real, Y(u> n ) is complex, and 
both are linear functions of B. For arbitrary complex 
gap functions A A (fc), Eqs. (|5"3"|) can be transformed into 
a system of nine equations for the components of X(uj n ), 
KeY(uj n ), and ImY(w„). 
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Plugging the expressions (|5Tj) and (152} into Eq. (08}, 
we obtain: 



f d?k 



n X 



d 3 k 

-{GxF* x )Y 3 + (G X F X )Y*] 
iTrV B N F TTj2[Xi(un)-X ,i(u}„)] (54) 



[here we also used the first of the equations ([5"3"}]. The 
vertex correction to the susceptibility is then found from 
X2,ij = dM^/dBj. 

The total spin susceptibility has the form Xij — Xi.ij + 
X2,ij = Xij+xtj+X7j+X2,ij- Putting together Eqs. (09} 
and (|54p . we arrive at the following final expression: 

EdXi (oj n ) 
dB — , (55) 

n J 

where XN.ij is the normal-state susceptibility (|27p . which 
is not sensitive to disorder. At T — > 0, we have X(uj n ) — > 
Xo(w„), with Cb n — uj n and D\ = A^. Using the identity 



|A A (fc)|= 



( W 2 + |A A (fe)P)3/2 



l-Y x (k,T), 



where Y x is given by Eq. (|23|) . we recover the suscepti- 
bility (|26|) of a clean superconductor. 



IV. RESIDUAL SUSCEPTIBILITY 

The equations derived in the previous section cannot 
be solved analytically in the general case. On the other 
hand, application of our results to real noncentrosymmct- 
ric materials is complicated by the lack of a definite infor- 
mation about the superconducting gap symmetry and the 
distribution of the pairing strength between the bands. 
Therefore one can only make progress by using some sim- 
ple models. 

First, we assume that the pairing corresponds to the 
unity representation of the point group and is fully 
isotropic: 4>+(k) = <p-(k) — 1, i.e. A x (k) = r)\. Then it 
follows from Eqs. (|46} that D x (k,uj n ) = D(oj n ) and the 
order parameter components (which can be chosen to be 
real) satisfy the self-consistency equations: 

EV^=T^J-^F x (k,» n ) 



■Kp X N F Ty^^ 



D\(uJ n ) 



(56) 



where ui n is found from Eq. (|4"S"j) . The Matsubara sum is 
cut off at the limiting frequency of the order of the BCS 
shell width e c . 

We further assume that the SO coupling is weak com- 
pared to the Fermi energy, |7|/ef "C 1, and that the 



pairing strength, see Eq. (fT0|) does not vary between 
the bands: V ++ = V__ > 0. Since the matrix of the 
coupling constants is symmetric and positive-definite, we 
also have U+_ = U_+ and V++ > \V+-\. The self- 
consistency equations have two solutions: 1) the order 
parameter magnitudes and phases in the two bands are 
the same: i] + = r/_ = rj, and 2) the magnitudes are the 
same, but the signs are opposite: ?y+ = — rj- = rj. Ac- 
cording to Sec. Ill Bt the former solution corresponds to 
the singlet state in the spin representation, while the lat- 
ter - to the "protected" triplet state. While in the clean 
limit the spin susceptibility for both states is given by 
Eq. (|29|) . the effects of impurities in the two cases have 
to be analyzed separately. 



A. 7]+ = 7]- = rj 

In this case D + = Z?_ = D(oj n ), and Eqs. (|4"5|) and 
46)) take the following form: 



D = f] + T- 



D 



l^n = W„ + r- 



The solution of these equations is D{tu n ) = Z{uo n )rj, ij n = 
Z(u> n )u> n , where Z{uj n ) = 1 + T/*Ju% + rj 2 . The gap 
equation ([56"} becomes 



n 9i 



D 



= , (57) 



where gi = (V++ + Vl\ )N F is the dimensionless coupling 

constant. Note that the scattering rate has dropped out 
of the gap equation, so that there is an analog of the 
Anderson theorem: neither the gap magnitude nor the 
critical temperature are affected by impurities: T c (r) — 
T c i = (2je c /n)e- 1 ^ 1 (In 7 ~ 0.577 is Euler's constant). 
In particular, the gap magnitude at T = is given by 
the clean BCS expression: r)(T = 0) = 770 = {i: / "i)T c \. 

Let us now calculate the spin susceptibility, using Eq. 
(155)) . Since the gap functions are real, Y is purely imagi- 
nary: Y(uj n ) = iY(uj n ), and Eqs. ([55} take the following 
form: 



1 - Ci 

C-2 1 — C's 



X 
Y 



HbC 2 B 

where C m ^ = C m (^j) F , m = 1,2,3, and 
Ci«)=r 

c 2 (w„) = r 



(58) 



c 3 (w n ) = r 



1 





Z{u n ) (ujI 


+ ^2)3/2 


1 




Z(uj n ) {wl 


+ ^2)3/2 


1 





Z( Un ) (^+^)3/2 
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calculated analytically, and we obtain the results plotted 
in Fig. M 

Thus we see that, similarly to spin-orbit impurities in 
a usual centrosymmetric superconductor^ scalar impu- 
rities in a noncentrosymmetric superconductor lead to an 
enhancement of the spin susceptibility at T = 0. Since 
the interband contribution is not sensitive to disorder, 
this effect can be attributed to an increase in the intra- 
band susceptibilities. 



B. 



v+ = -v- = n 



17 770 



FIG. 2: The residual susceptibility vs disorder strength for 
77+ = r?_ = 770. The solid line corresponds to the trans- 
verse components in the 2D case {Xzz = Xp an d i s disorder- 
independent), the dashed line - to all three diagonal compo- 
nents in the 3D case. 



Substituting the solution of the equations (|58|) in Eq. 
(|55|) . we obtain: 



Xij 



= XP$ij - Kfi 2 B N F T^2 A y ( 59 ) 



where A = 4r(l - L) _1 L, and L = 61 + 62(1 - C 3 y 1 C 2 . 
In the coordinate system in which (7i7j)_F is diagonal, 
the susceptibility tensor is also diagonal: 



Xu(T) 
XP 



(60) 

where I\ = (1 - (j?)f)T. 

We are particularly interested in the effect of disorder 
on the residual susceptibility at T = 0. In this limit, the 
Matsubara sum in the last formula can be replaced by a 
frequency integral, which gives 

MT=0) = 1 m F + (*) , (6i) 

Xp Vw 



where 



* l( *) = l-_ 1- 



: arctan 



1-x 
1 + x 



(62) 



[at x > 1 this function is evaluated using arctan(ix) = 
i tanh _1 (;c)]. While the first two terms on the right-hand 
side of the expression (|6ip represent the residual suscep- 
tibility in the clean case, see Eq. (|28|) . the last term de- 
scribes the impurity effect. In a weakly-disordered super- 
conductor, using the asymptotics < &i(x) ~ 7rx/4, we find 
that the residual susceptibility increases linearly with dis- 
order. In the dirty limit, T 3> r]o, we have <&i(x) — ► 1, 
therefore Xa(T = 0) approaches the normal-state value 
Xp- For the two simple band-structure models discussed 
in the end of Sec. Ill Al the Fermi-surface averages can be 



In this case D + = —D- = r), and we obtain from Eqs. 
451 and (l56l): 



W„ = UJ n + r 



(63) 
(64) 



where <72 = )Np. In the absence of impurities, 

the critical temperature is given by the BCS expression: 
T c (r = 0) = T c2 = (2 7 e c /7r)e- 1 / 92 . If > (attrac- 
tive interband interaction), then g 2 < g\ and T c2 < T c \, 
i.e. the phase transition occurs into the state 77+ = ry_ . If 

V-\ < (repulsive interband interaction), then g 2 > gi 

and T c2 > T c \, i.e. the phase transition occurs into the 
state 77+ = — rj-. 

In contrast to the previous case, both the critical tem- 
perature and the gap magnitude are now suppressed by 
disorder. Indeed, linearizing the equations (j6"3"|) and (j6"4"|) 
with respect to 77, we obtain the following equation for 
the critical temperature T c (r): 



2irT r J \ 2 



(65) 



Thus the suppression of T c by scalar impurities is de- 
scribed by the same Abrikosov-Gor'kov function as in 
a conventional BCS superconductor with magnetic im- 
purities. The superconductivity is completely destroyed 
if the disorder strength exceeds the critical value r c = 
(n/2 7 )T c2 . 

To find the gap magnitude at T = we follow the 
procedure outlined in Ref. [jOl Replacing the Matsubara 
sum by a frequency integral in the gap equation (|64[) . we 
have 

l-fduY' —rL-) 

92 Jo \ \JCj 2 + rj 1 y/uj 2 + ff I 



duj- 



1 



!! \/u 2 +jf 



(66) 



where the function Q(uS) is obtained from the continuous 
limit of Eq. 



Q = uj + T- 



y/uj 2 + T] 2 



(67) 
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Since e c ^> r), the last integral in Eq. (|66|) is equal to 
ln(2e c /rf). In the clean case we recover the BCS expres- 
sion for the gap magnitude at T = 0: 770 = {^/"f)T C 2 = 
2T C . In the presence of disorder Eq. (j6"6"|) can be repre- 
sented in the form 



In* 



1 



1 



\JtjJ 2 + rj 2 y/Cj 2 + T] 2 



(68) 



Using Eq. (|67|) the second term can be transformed into 
an integral over Co. Introducing the notation x = T/r), we 
arrive at the following equation for the gap magnitude at 
zero temperature as a function of T: 



(69) 



where 



7TX 

T{x) = — - hx{2x) + 6(x - 1) 



ln(a; + \J x 2 - 1) 



- — arctan \J x 2 — 1 — 
2 



2x 



The equation 1|69|) does not have solutions at T > T c , 
which is consistent with the complete suppression of su- 
perconductivity above the critical disorder strength. 

Repeating the steps from the previous subsection, we 
obtain the susceptibility in the form (JHH), where 



Lij(uJ n ) — T— 



if 



{Cb 2 n +7,2)3/: 



■(7»7i) 



In the coordinate system in which (7i7j)i? is diagonal, 
the nonzero components of the susceptibility tensor are 
given by 



Xu{T) 
Xp 



l-( 7 f) F7r T£- 



(70) 

We note that for a spherical 3D model with (7?) p = 1/3 
this expression has exactly the same form as the suscep- 
tibility of the superfluid 3 He-B in aerogel, see Refs. l3ll 
and HE 

At zero temperature we replace the Matsubara sum 
by a frequency integral and then transform it into an 
integral over Co using Eq. (|67[) . with the following result: 



Xu(T = 0) 
Xp 



(% 2 



$ 2 _ 



1 F 'i>2 



(71) 



where 



* 2 (ar) 



(y 2 
1 



I) 3 / 2 



(y2 + 1)3/2 _ x{ A f 2 )l 



(72) 



and y m j n = 6{x — \ )\/ x 2 — 1. The last term on the right- 
hand side of Eq. J7T]) describes the effect of impurities 




r/)?o 

FIG. 3: The residual susceptibility vs disorder strength for 
?7+ = — t?_ (r/770 = 0.5 corresponds to the critical disorder 
strength V — T c ). The solid line corresponds to the trans- 
verse components in the 2D case (xzz = Xp an d i s disorder- 
independent), the dashed line - to all three diagonal compo- 
nents in the 3D case. 



on the residual susceptibility. For each disorder strength, 
one first finds the solution x(T) of Eq. (f6"9"|) and then cal- 
culates <i>2(^(r)). In the weak disorder limit we have 
x{T) ~ T/2r c « 1, and $ 2 (z) ~ (37re/16)(l - ( 7 2 )f), 
i.e. the residual susceptibility increases linearly with dis- 
order. At T -> T c = ryo/2: x(T) ~ ^r c /12(r c - f) > 1. 
In this limit $ 2 (ie) — > 1 and Xm( t = 0) — > Xp- 

The dependence of x*i(T = 0) on the disorder strength 
for the 2D and 3D models is plotted in Fig. 03 As in the 
case 77+ = 77_, the residual susceptibility is enhanced by 
impurities. 



V. CONCLUSIONS 

The spin susceptibility of a clean noncentrosymmet- 
ric superconductor has a large residual value at T = 0, 
which can be attributed to the temperature-independent 
contribution of the virtual transitions between the non- 
degenerate bands split by the SO coupling. We have 
studied the effects of disorder and found that scalar im- 
purities in noncentrosymmetric superconductors act like 
spin-orbit impurities in centrosymmetric BCS supercon- 
ductors, in the sense that they considerably enhance the 
residual susceptibility. The quantitative details depend 
on many factors, in particular the shape of the Fermi 
surface, the structure of the SO coupling, characterized 
by 7(fe), and the symmetry of the order parameter. 

In the model in which the gap functions in both SO 
split bands are the same and isotropic, the critical tem- 
perature T c does not depend on the elastic scattering 
rate T, while the residual susceptibility increases and ap- 
proaches the normal-state value xp m the dirty limit 
r — ► 00. If the two gap functions have the same mag- 
nitudes but opposite signs, T c is completely suppressed 
above the critical concentration of impurities correspond- 
ing to r c , with the residual susceptibility approaching \p 
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at r -> r c . 

The band structure models considered here cannot be 
directly applied to CePtaSi, whose Fermi surface is quite 
complicated and consists of multiple sheets^ It is not 
known which one (or ones) of them are superconduct- 
ing. The order parameter symmetry is not known ei- 
ther, although there is strong experimental evidence that 
the superconducting order parameter has lines of gap 
node o 33 ' 34 i 35 ' 36 i 37 and likely corresponds to one of the 
nontrivial one-dimensional representations of the point 
group C4,,. We expect however that the qualitative pic- 
ture described above remains valid for anisotropic pairing 
as well, namely the residual susceptibility increases in the 
presence of impurities. 
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APPENDIX A: SPIN-ORBIT COUPLING OF 
BAND ELECTRONS 

In order to derive the expression ([T]), we start with 
the Hamiltonian for non-interacting electrons in a perfect 
crystal lattice: 



n 2 fr 2 



2m 



U(r) 



Am 2 c 



:&[VU(r) x V], (Al) 



where U(r) is the lattice potential. The last term, Hso, 
represents the spin-orbit coupling. In the absence of the 
latter, the eigenstates of H are the Bloch spinors: 

(ra\kpa) = — w. klJi (r)e ihr x a {a), (A2) 
v V 



where V is the system volume, a is the spin projection, 
tyhn^y) nave the same periodicity as the crystal lattice, 
and Xa are the basis spinors: Xa( tJ ) — The Bloch 
spinors are labelled by the wave vector k, the band index 
/j, and the spin index a, which distinguishes the states 
within the same band. The eigenvalues e At (fc) describe 
the electron dispersion in the /xth band and have the 
following symmetry properties: e^(fe) = e (tt (— fe), = 
e^(g k), where g is any operation from the point group 
of the crystal. 

Next, we calculate the matrix elements of H$o hi the 
basis of the Bloch states (IA2I): 



(k^a\Hso\pvP) 



£ 

ijk 



J d 3 r6 ]k (r)e^P-V r 7 (A3) 



where i, j, k = x,y, z, and 

e jk = {VjU)<p* kfJt (-iVk+Pk)<ppv 

Since Qjk are lattice-periodic functions of r, the integral 



in Eq. (|A3[) is nonzero only if p — k = G, where G is 
a reciprocal lattice vector. Because both k and p are in 
the first Brillouin zone, the only possibility is p = k. The 
Hamiltonian remains a nondiagonal matrix in both the 
band and spin spaces, and has the following form in the 
second-quantization representation: 



Ha 



k,fiv,a.{3 



^(tyS^Sap + r y ilv {k)(T a p\al fJtot a kv0 , (A4) 



where the chemical potential is included in the band dis- 
persion functions, and 



Am 2 c 



2„2 



x- [ d 3 r [(VU) x ^ H V + k)<p kl/ ] (A5) 

^ J V 

(the integration is performed over the unit cell of volume 
v). The expression (|A4[) is exact for non-interacting elec- 
trons, regardless of the band structure and the strength 
of the SO coupling. 

The functions 7u„(fe) satisfy certain symmetry- 
imposed conditions. Since Hq is Hermitian, we have 



(A6) 



Under the point group operations g, e.g. rotations, 
the second-quantization operators transform as follows: 
<*U 7* Epat^fPM where U{g) = e -^™)l 2 is 
the spinor representation of the rotation about a direc- 
tion n by an angle 9^ Requiring that Hq remains in- 
variant under g we obtain: 



(A7) 



Under time reversal, fa' h/ia -> /* X^( icr 2)a/3al fcAl/3 (/ is 
an arbitrary c- number coefficient), therefore "^(fe) — > 
— 7* ly (— fc). If the time-reversal symmetry is not broken, 
we have 



l^{k) = -% v {~ k )- 



(A8) 

a , and 7 M „ trans- 



Finally, under inversion a k - 
form like pseudovectors: 7 A1 „(fc) — > 7 M „(— fc). 

In a centrosymmctric crystal, 7 /1I /(fc) = 7 M1 ,(— k), 
therefore, using Eq. (|A8[) . 

lUk) = ~%u(k). (A9) 

It follows from the conditions (|A6|) and (|A9|) that 
l^^{k) = 0. Therefore one needs to include at least two 
bands in Eq. (|A4|) . in which case 712 (k) = — 721 (fe) = 
i£(k), and the Hamiltonian takes the form 

Ho = € ^ k )al lia a k ^ a 

k,a fi=l,2 

+i £(k)(T a0 (ai la a k2f3 - a{ 2a a kll3 ), (A10) 

k,a0 
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see also Ref. H^. The pseudovector I is real, even in fc, 
and satisfies l{k) = g£(g" 1 fc). 

In contrast, in a noncentrosymmetric crystal, the con- 
straint (|A9|) is absent, and the effects of SO coupling can 
be studied in a minimal model in which one keeps just 
one band in the Hamiltonian (|A4j) . The band index fj, 
can then be dropped and the SO coupling can be de- 



scribed by a single pseudovector function 7(fc), which is 
real, odd in fc, and invariant with respect to the point 
group operations: 7(fc) = g , y{g^ 1 k). In this way one 
arrives at the effective band Hamiltonian |T]). Although 
in principle one can calculate 7(fc) using Eq. (|A5[) . we 
consider it as a model parameter. 
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